Two-Photon Scattering in One Dimension by Localized Two-Level System 
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We study two-photon scattering in a one-dimensional coupled resonator arrays (CRA) by a two-level system 
(TLS), which is localized as a quantum controller. The S -matrix is analytically calculated for various two- 
photon scattering processes by TLS, e.g., one photon is confined by TLS to form a bound state while the other 
is in the scattering state. It is discovered from the poles of the S -matrix that there exist two kinds of three-body 
bound states for describing two bound photons localized around TLS. 
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Introduction. — For the architecture of a new generation all- 
optical quantum devices, we need to study various physical 
mechanisms of the coherent photon transports in the low di- 
mensional confined structure with sub- wavelength scale. Re- 
cently, the single photon transmissions in a waveguide and one 
dimensional (ID) coupled resonator arrays (CRA) with some 
controllers, e.g., a two-level system (TLS), have been exten- 
sively studied iil,l2l[3l|4|,[5|]. It * s worthy to notice that, in such 
CRA hybrid system, there usually exist bound states of single 
photon around the localized TLS , thus TLS can behave 
as a quantum controller to realize a quantum switch [6] for 
photon transports. 

We emphasize that most of investigations in this area fo- 
cused on the cases with single photon, and two-photon case 
was only considered for the linear dispersion relation of 
waveguides In this Letter, we develop a two-photon 

scattering approach for the general case with arbitrarily given 
nonlinear dispersion relation. Actually, the effects of TLS on 
the quantum natures of multi-photon statistics, such as photon 
bunching and anti-bunching, have not been investigated sys- 
tematically for the nonlinear dispersion relations. However, 
a comprehensive understanding for fundamental processes of 
multi-photon scattering is necessary to both practical applica- 
tions and theoretical explorations. In experiments, the prac- 
tical processes may concern two or more photons, which ob- 
viously affect on the efficiency of single photon creation and 
transfer. In theoretical studies, the hybrid system with many 
photons is related to the Lee model [80 in the sector with high 
excitations. 

Our hybrid system consists of the CRA system (Fig. Hk) 
with a nonlinear dispersion relation (3, fl 0] and a TLS local- 
ized in one of the resonators. It has been shown that in some 
circumstances there exist two single-photon bound states, and 
energy band structure for a scattered photon, which are dis- 
played in Fig.[T[). With this illustration, we address our main 
results: (a) when scattered by TLS, which simultaneously 
binds another photon to form a bound state, the incident pho- 
ton is elastically scattered under a condition we will discuss 
as follows; (b) when two incident photons are scattered by 
TLS, the photon correlation is induced, and the nonlinear dis- 
persion results in richer quantum statistical characters of out- 
going two photons, which can be controlled by TLS; (c) there 
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FIG. 1 : (Color online) (a) The schematic for the CRA system: the red 
circle denotes the two level system (TLS), and the blue dots denote 
the coupled resonators, (b) The band structure of the CRA system for 
the single excitation case: the continuum displays the single photon 
scattering states while the blue lines above the band top and below 
the band bottom represent the single photon bound states. 



exist two kinds of three-body bound states for describing two 
bound photons localized around TLS. 

Model — The model Hamiltonian of our hybrid system 
reads 

H = Q \e)(e\ + J dks k a\a k + V J -^{a\cr- + H.c.) , (1) 

where the operator cr~ = \g)(e\ denotes the flip from the 
ground state (GS) \g) to the excited state \e) with the energy 
level spacing Q. The nonlinear dispersion relation of photons 
in CRA is s k = ojq - 2 J cos k, where ojo is the eigen-frequency 
of each cavity, / is the hopping constant characterizing the 
inter-cavity coupling in the tight-binding approximation, k 
is the momentum of photon, and the inter-cavity distance is 
taken as unity. Here, a k (a|) is the annihilation (creation) op- 
erator for the photon with momentum k, and V is the coupling 
constant of TLS and photon. 

In the single excitation subspace spanned by the basis 
{|0) \e) , al |0) \g)} JES], the system possesses an energy band 
of width 4 J centered in ojq and two single-photon bound states 
\E S ) = V^[|0>k> + V(27T)-^ 2 f dk(E s - e)- l al\0)\g)] (s = 
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FIG. 2: (Color online) (a) The schematic for the photon scattering 
by TLS in the single photon bound state: the first line indicates the 
single photon incident, the second line means the scattering into the 
out-going state of single photon, and the third line shows that the 
bound state is broken so that there are two out-going photons; (b) 
the schematic for two-photon scattering by the TLS in GS: the first 
line indicates the two photon incident, the second line means the 
scattering into two out-going photons, and the third line shows that, 
while a bound state form, another photon is scattered into out-going 
state. 



1, 2) with energies E\ below the band bottom s m [ n = cjq - 2 J 
and E2 above the band top £ max = ojo + 2J (Fig. \^p). Here, 
the normalization constant is Z" 1 = 1 + V 2 J dk(E s - Sk)~ 2 /2n 
(s = 1,2). E s are two real solutions of G~ l (E s ) = HH, 
where G(oj) is the Green's function of photon, and G~ 1 (oj) = 
co-Q- 1,(0)). Here, Y(oj) = V 2 f dk(co - £ k )~ l /2n is the self 
energy of photon. 

For two-photon scattering, there may exist three possible 
processes: (a) one incident photon with momentum k scat- 
tered by TLS in the bound state \E\) or |#2) (see Fig.[Tk); (b) 
two incident photons scattered by TLS in GS; (c) two photons 
both bound by TLS form the three-body bound state. We point 
out that the behaviors of two photons in CRA are essentially 
different from those of photons propagating in a waveguide, 
in which the processes (a) and (c) do not exist due to its linear 
dispersion relation. The processes (a) and (b) are schemati- 
cally shown in Fig.|2^ andEJ), respectively. 

Scattering eigen-states. — The scattering eigen-states for the 
two-photon processes are assumed to be 



|¥ (/) } = J dk®f(k)a\ |0> \e) + J dkdk?®®(k,k?)ala 

(2) 

for / = 1,2,3. The corresponding eigenvalues are = 
Ei + £k for / = 1 , 2 and £ (3) = + Sk 2 • Here, the eigen- 
states |*F (1 ' 2 ^ describe one photon of momentum ko scattered 

by TLS in the bound state ^1,2), while |^ (3) ^ describes the 
two photons with momenta k\ and ki scattered by TLS in GS. 
For convenient, we set so = Sk Q , £i = £ = £k and s' - £# 
below. 

According to the three-particle scattering approach 
llUl2h for the Lee model [8D, the solutions of scattering eigen- 



states for the secular equation H \^ s) ) = E (s) |Y< 5 >) (5=1,2) 
are obtained as 



V[Og(Q + 0< J) (*Q] 
2V2tt(£ ( * ) -£-£ / + /0 + )' 



(3) 



where & e s \k) = ^/Z~ s [6kk + i/f(s,so)] contains the incident 
component Skk Q in the k- space. The scattering part for single 
photon reads as 



iff(e, so) 



V 2 



£0 — Eo 

[^-^G + (e ) 



27T(£o — e) e - E 

2A s (£,£q)G + (£q) 
Z s G + (£o)-A s (£ ,£ ) 



(4) 



where A s (e,e ) = Z- S (E- S - E S )G + (E& - E,)/(E^ - E- s - 
s) + f max ds f J + (E^ - s\ s)(s' - E s )ImG + (s'), and we define 
F+(x) = F(x + /0 + ) for arbitrarily given function F(x). Here, 
J(E, e) = 7T- l G(E)/(s - E) and s ± s (s = 1,2). 

For two-photon scattering by TLS in GS, the two-photon 
wave-function 



& 3 \k,k') = dk kl Sk>k 2 + 



V[<S>?\k) + O^Q] 
2 V2^(£ (3) -s-s' + /0 + ) 



(5) 



contains the incident two-photon component 5kk l 8k f k 2 - an d the 
two-photon wave-functions in the scattering part 
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[V 2 (P(£)G + (£i)G + (£ 2 )+2nG + (E^-£) ^ S kkl l 

i 

(6) 

where 0(e) = 21(e) /A - £ !=u (e - £;)"'• Here, 1(e) = 
Z s=h2 Z s G + (E< 3 > - E,)KE® -E s -e)+ f £m ~ de'J + (E^ - 
£',e)ImG + (e'), A = C - >Z S=U2 Z S G+(E (3) - E s ), and C = 
tC x f n,ax de'G + (E (3) - e')ImG + (e'). 

Scattering by TLS in bound states. — For an incident photon 
with momentum ko scattered by TLS in the bound state \E S ), 
the out-going state |out> contains the single photon component 
a ?Q \E S ) and the two-photon out-going part a Pl a ?2 10) \g). The 
above obtained scattering eigen-states |^ L2 ^ result in the S- 
matrix iHHEllll. Its element 
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represents an elastic process "7 + BS —> y + BS ": a photon 
y with momentum &o scattered by TLS in the bound state \E S ) 
into a photon y with momentum po while the bound state is 
unchanged. The transmission and reflection coefficients are 



t[ s) = 1 + r [ s) and 

k Q k 



(s) _ iV 2 G + (£o) Z s G + (£ )+A s (£ ,£ ) 
r ^ 2J |sin M Z s G + (£ ) - A s (£ , e ) ' 

Another S -matrix element 



(8) 
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PiP2ko = i2Z s G + (£ )G + (£ Pl )G + (£ P2 ) 
6(E^-£ P1 -£ P2 ) 



Z s G + (£ )-A s (£ ,£ ) 



(9) 



3 




(c) (d) 




(c) (d) 



FIG. 3: (Color online) The transmission and reflection probabilities 
and outgoing wave-function of two photons: the parameters of the 
system are ojq = 4 and J - V = 1. In (a) and (c) Q = 2.5, and 
in (b) and (d) Q = 5. The reflection probability, (blue dotted line) 
transmission probability (red dashed line) and the sum (green solid 
line) of them are drawn in (a) and (b) for the photon scattered by TLS 
in the lower and higher bound states, respectively. The coordinate 
representation of two photon out-going wave-functions are shown in 
(c) and (d) for photon with energy so = 4.65 scattered by TLS in 
the lower bound state and photon with energy s = 2.55 scattered by 
TLS in the higher bound state, respectively. 



means an inelastic process "y + BS — > 2y": a photon y with 
momentum & scattered by TLS in the bound state \E S ) into 
two photons with momenta p\ and p2. 

Let us analyze the physical processes described by Eqs. © 
and ©. When e < = 2s m \ n - E\, ^ vanishes, 
thus the incident photon is elastically scattered by TLS in the 
bound state \E\). When £q > s^, the bound state \E\) is 
stimulated by the incident photon, and the out-going two pho- 
tons emerge simultaneously. The incident photon with energy 
£ ° > £ ± ~ 2 £ max - E2 is elastically scattered by TLS in the 
bound state IE2) while the bound state IE2) is not affected. 
When so < s^? the out-going two photons emerge. 

As shown in Fig. [3^ and Fig.[3j), the reflective probabilities 

I n |2 

associated to the bound states are suppressed in com- 
parison with that for the single photon scattering in Ref. 01. 

Here, the probabilities pk = |r^| + |^| are also displayed 

in Fig. [3^ for s = 1 and Fig. [3]3 for s = 2: when so < 

(2) 

(so > s th ), the probability pk = 1 means the elastic scat- 
tering by TLS in the bound state \E\) QE2))', when sq > 
(so < Sfa), the probability pk Q < 1 means the inelastic scat- 
tering by TLS in the bound state \E\) QE2)), and in this case 
the two photons emit. The probabilities \Q> s (xi, X2)\ 2 for two- 



FIG. 4: (Color online) The second order correlation functions: the 
parameters of system are £1 = ojq = 4, V = 1 and ./ = 1, (a) the 
second order correlation function for the two photons both in reso- 
nance with TLS; (b) the second order correlation function for e\ = 3 
and s 2 = 5; (c) for s\ = s 2 = n, the dashed (Green) and solid (red) 
curves denote the second order correlation function when x c =0 and 
x c -l, respectively; (d) for E\ = 3 and s 2 = 5, the dashed (Green) and 
solid (red) curves denote the second order correlation function when 
x c -Q and x c =1.6, respectively. 

photon creation in the out-put, which is defined by 

0,(xi,jc 2 ) = J Sf x ^{^ XlM ^ + e ip2Xl+ipiX2 ), 

(10) 

are plotted in Fig. |3t and Fig.[3]i The numerical calculations 
show that the out-going two photons prefer to occupy around 
TLS. 

Two-photon scattering by one TLS. — For the two inci- 
dent photons both not bound by TLS, the out-going state 
|out) generally contains the single-photon part |one> = 
I* P , s =i2Spik 2 a l\ E s) and the two-photon part |two) = 
I^PiP2 S PiP2kik 2 al l al 2 10> \g). The Lippmann-Schwinger ap- 
proach for l^ 3 ^ gives the S -matrix explicitly, which elements 
are sj^ = -iV 3 W(2nZ s )-^ 2 S(E s + s p - £ (3) ) (s = 1,2) 
for the inelastic process "2y — > y + BS " and S PlP2 kik 2 = 
Sin + RS(E^ - s Pl - £ P2 ) for the elastic process "2y — » 2y". 
Here, W = G + (s p )G + (si)G+(s 2 )/A, and the first term S in = 
S Plkl S P2k2 + S P2kl S Pl k 2 of S PlP2kl k 2 describes the factorization 
of the two-photon scattering, where S p k = tkS(p-k)+rkS(p+k) 
contains the transmission coefficient tk = 1 + rfc and reflection 
one 



2J\sink\(s k -Q) + iV 2 ' 

They describe the single photon scattering for TLS in GS 01. 
The correlated photon scattering induced by TLS is described 



4 



by the second term of S PlP2 k v 



estimated as B\ 2.82 and B2 - 12.54, respectively. Here, 



(12) 



i=l,2 



which the exhibits the background fluorescence of two pho- 
tons AS 01. 

In the inelastic scattering of two photons, a part of incident 
photon energy is absorbed by TLS, and thus the out-going 
state |one) describes that one photon forms a bound state while 
another emits. Another out-going state |two) means that the 
two photons are elastically scattered by TLS without any en- 
ergy loss. When E\ + £ max < £ (3) < £2 + £min the inelastic 
process is forbidden and only the elastic scattering takes place. 
The second order correlation function [4] for two photons is 
G(xi,x 2 ) = \g(xi,x 2 )\ 2 , where 



g(xi,x 2 ) 



f 



dp\dp 2 
4n 



S pip 2 hk 2 ( e 



,ipixi+ip 2 x 2 _j_ e ip 2 x\+ip\x 2 



2 ) (13) 



is the out-going wave-function in the coordinate representa- 
tion. 

Two correlation functions G(x\, x 2 ) are drawn in Fig.|4^ for 
s\ = £ 2 = Q and in Fig.|4b for s\ = 3 and £2 = 5. In fact, due 
to the nonlinear dispersion relation, the correlation function 
G(x\,X2) not only depends on the center-of-mass coordinate 
x c = fe+*i)/2, but also on the relative coordinate x r = x<i~x\ 
of two photons. For the different x c the correlation functions 
G are shown in Fig. |4j: for s\ = £2 = £1 and in Fig. [4jl for 
£\ = 3 and £2 = 5. These figures indeed show the different 
photon statistics for the different center-of-mass positions of 
photons. 

Three-body bound states. — It follows from the poles of S- 
matrix that there exist two three-body bound states with en- 
ergies B s - E s + J3 S (s = 1,2), where J3 S are determined by 
Z S G+(J3 S ) = A S (J3 S ,J3 S ). It can be proved that J3\ < £ m { n and 
Pi > £max- Using the approach in Ref. [15], we obtain the 
bound state 

\B S ) = N[ J dkrfi\k)a\ |0>k>+ J dkdk f Tf s \k,k f )a\al |0> \g)] 

(14) 

(s) 

by the secular equation H \B S ) = B s \B S ). Here, rf e (k) = 
V 2 (J3 S - sT l A s (s,B s )/27T, 



2V2ti(B s -£-£') 



(15) 



and N is the normalization constant. The wave-functions 
n ( e s \x) = f dkri { e\k) exp(to)/ ^^?JI and 



// !l (. f| ,.v 2 ) - I ^^rf s) {k u k 2 ){e ipiXl+ip ^- + j™*™) 

(16) 

in the coordinate representation both tend to zero, when x, 
x\ and X2 tend to infinity, respectively. For the parameters 
a>o = 4, 7=1, Q = 3, and V = 2, the binding energies are 
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FIG. 5: (Color online) The wave-functions |^(x)| 2 for s = 1,2: the 
parameters of system are coq = 4, J = 1, Q. = 3, and V = 2. As the 
photon goes to infinity, the probabilities tend to zero. 



|^(x)| are shown in Fig. [5] for s 



1,2, which both go to 
zero as x — > oo. 

Conclusion. — We have investigated the two-photon scatter- 
ing in the hybrid CRA system. The two-photon S -matrix are 
calculated for various scattering processes by TLS in the CRA 
system with nonlinear dispersion relation. We show that the 
three-body bound states exist for describing the two bound 
photons localized around TLS. The binding energies are de- 
termined explicitly by the poles of S -matrix. For the two- 
photon scattering by TLS in GS, the quantum statistics of 
scattered photon are analyzed in detail. The similar approach 
can also be applied to dealing with scattering problems of the 
TLS-photon hybrid system with more complicated architec- 
tures and dispersion relations. 

One (T. Shi) of the authors would like to thank Q. Ai 
for many helpful discussion. The work is supported by Na- 
tional Natural Science Foundation of China and the National 
Fundamental Research Programs of China under Grant No. 
10874091 and No. 2006CB921205. 



[1] J. T. Shen and S. Fan, Phys. Rev. A 76, 062709 (2007). 

[2] L. Zhou, Z. R. Gong, Y. X. Liu, C. P. Sun, and F. Nori, Phys. 

Rev. Lett. 101, 100501 (2008). 
[3] D. Z. Xu, H. Ian, T. Shi, H. Dong, and C. P Sun, 

larXiv:0812.0429l 
[4] T. Shi and C. P. Sun, Phys. Rev. B 79, 205111 (2009). 
[5] J. T. Shen and S. Fan, Phys. Rev. A 79, 023837 (2009). 
[6] C. P Sun, L. F. Wei, Yu-Xi Liu, and F. Nori, Phys. Rev. A 73, 

022318 (2006). 

[7] J. T. Shen and S. Fan, Phys. Rev. Lett. 98, 153003 (2007). 
[8] T. D. Lee, Phys. Rev. 95, 1329 (1954). 

[9] M. L. Goldberger and S. B. Treiman, Phys. Rev. 113, 1663 
(1959). 

[10] C. M. Sommerfield, J. Math. Phys. 6, 1170 (1965). 

[11] R. D. Amado, Phys. Rev. 122, 697 (1961). 

[12] A. Pagnamenta, J. Math. Phys. 6, 955 (1965). 

[13] G. C. Wick, Rev. Mod. Phys. 27, 339 (1955). 

[14] G. Kallen and W Pauli, Dan. Mat. - Fys. Medd. 30, 7 (1955). 

[15] A. Pagnamenta, J. Math. Phys. 7, 356 (1966). 



